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Abstract 

We examine various generalizations e.g, exactly solvable, quasi exactly solvable and non Hermitian variants 
of the quantum nonlinear oscillator. For all these cases the same mass function has been used and it has also 
been shown that the new exactly solvable potentials posses shape invariance symmetry. The solutions are 
obtained in terms of classical orthogonal polynomials. 



1 Introduction 

Recently there has been a surge of interest in obtaining exact p] and quasi exact solution [2] of the position 
dependent mass Schrodinger equation (PDMSE) for various potentials and mass functions by using various 
methods like Lie algebraic techniques [5] , supersymmetric quantum mechanics (factorization method) [H [5] , 
shape invariance approach [6] , point canonical transformation [7j , path integral formalism 8 , transfer matrix 
method [9] etc. Apart from the intrinsic interest the motivation behind this issue arises because of the relevance 
of position dependent mass in describing the physics of many microstructures of current interest, such as 
compositionally graded crystals [TU], quantum dots [TT], 3 He clusters [12], metal clusters [12] etc. The concept 
of position dependent mass comes from the effective mass approximation [14) which is an useful tool for studying 
the motion of carrier electrons in pure crystals and also for the virtual-crystal approximation in the treatment 
of homogeneous alloys (where the actual potential is approximated by a periodic potential) as well as in graded 
mixed semiconductors (where the potential is not periodic). The attention to the effective mass approach stems 
from the extraordinary development in crystallographic growth techniques which allow the production of non 
uniform semiconductor specimen with abrupt heterojunctions. In these mesoscopic materials, the effective 
mass of the charge carriers are position dependent. Consequently the study of the effective mass Schrodinger 
equation becomes relevant for deeper understanding of the non-trivial quantum effects observed on these 
nanostructures. The position dependent (effective) mass is also used in the construction of pseudo-potentials 
which have a significant computational advantage in quantum Monte Carlo method [15] . It has also been 
found that such equations appear in very different areas. For example, it has been shown that constant mass 
Schrodinger equation in curved space and those based on deformed commutation relations can be interpreted 
in terms of PDMSE in flat space [TB] and 'PT-symmetric cubic anharmonic oscillator [T7] . 
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The nonlinear differential equation 



(1 + Xx 2 )x - (Xx)x 2 + a 2 x = 0, A>0 (1) 

was studied by Mathews and Lakshmanan in [181 119) as an example of a non-linear oscillator and it was shown 
that the solution of ([I]) is 

x = Asin(ujt + 4>) (2) 
with the following additional restriction linking frequency and amplitude 

" 2 = TT\¥ < 3 » 

Furthermore (TTJ) can be obtained from the Lagrangian [TS] 

so that both the kinetic and the potential term depend on the same parameter A. So this nonlinear oscillator 
must be considered as a particular case of a system with a position dependent effective mass. Recently in 
a series of papers |20[ 12 1 j this particular nonlinear system has been generalized to the higher dimensions 
and various properties of this system have been studied. The classical Hamiltonian corresponding to the 
A-dependent oscillator is given by [TH1 [H] 



H =^) P ' + G) 5 (lWj' r* = VTTA^ Px , 9 = ma> (5) 

p x being the canonically conjugate momentum defined by p x = L being the Lagrangian and m is the 
mass. 

It has been shown in [5T] that in the space £ 2 (5ft, dfi) where dfi = ( , 1 ^^ i )dx, the differential operator 



„2 



y/1 + Xx 2 -^ is skew self adjoint. Therefore, contrary to the naive expectation of ordering ambiguities, the 
transition from the classical system to the quantum one is given by defining the momentum operator 

d 



P x = -iVl + Xx 2 — (6) 
dx 



so that 



^+^1 ) 



(l + Aa^K 

Therefore the quantum version of the Hamiltonian with h = 1 becomes |21j 

^=-^(l + Ax 2 )^- (±) Xx-^ + ~g(—^—) (7) 

2m v 'dx 2 \2m J dx 2* \l + Xx 2 J y ' 

where g = a(ma + A). It is to be noted that in ref [5T] the value of the parameter g has been slightly modified 
from that given in equation ([5]). 

It may be pointed out that this A-dependent system can be considered as a deformation of the standard 
harmonic oscillator in the sense that for A — * all the characteristics of the linear oscillator are recovered. 

In ref [21], the PDMSE corresponding to this nonlinear oscillator has been solved exactly as a Sturm- 
Liouville problem and A-dependent eigenvalues and eigenfunctions were obtained for both A > and A < 0. 
The A-dependent wave functions were shown to be related to a family of orthogonal polynomials that can 
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be considered as A-deformations of the standard Hermite polynomials. Also the Schrodinger factorization 
formalism, intertwining method and shape invariance approach were discussed with reference to this partic- 
ular quantum Hamiltonian. The existence of a A-dependent Rodrigues formula, a generating function and 
A-dependent recursion relations were obtained. 

In this paper our objective is to re-examine this problem and obtain closed form expression for the normalisa- 
tion constant, modified generating function and recursion relations for A(= ^)-deformed Hermite polynomials. 
A relation between the A deformed Hermite polynomials and Jacobi polynomials will also be obtained. We 
shall also obtain a number of exactly solvable, quasi exactly solvable and non Hermitian potentials corre- 
sponding to the same mass function m(x) = (1 + Ax 2 ) -1 . It will be seen that some of these potentials are 
generalizations of the nonlinear oscillator potentials while the others are of different types. It will be shown 
that these exactly solvable potentials are shape invariant. Moreover these potentials can also be complexi- 
fied and by doing so we shall also obtain a number of exactly solvable non Hermitian potentials within the 
framework of PDMSE. As a method of obtaining these results we shall use point canonical transformation 
consisting of change of coordinate only. The organization of the paper is as follows: in section[2]we shall obtain 
exactly solvable potentials and a relation between A-deformed Hermite polynomials and Jacobi polynomials; 
in section [3] it is shown that the exactly solvable potentials are shape invariant; in section [5] we obtain exactly 
solvable non Hermitian potentials; section [5] deals with complex quasi exactly solvable potentials and finally 
section [6] is devoted to a discussion. 



2 Exactly solvable potentials for the mass m(x) = ( i + \ x 2 ) 

Here we shall obtain exact solutions PDMSE for a number of potentials with the same mass function m(x) = 
1 ' For this purpose we first write the PDMSE corresponding to the Hamiltonian given in Eqn. ([7]) with 

dip g ( 1 



(l+Ax a ) 



m = 1 and A > as 



Xx 



dx A V 1 + Ax 2 



%p = Ei> 



E = 2e-\ 
A 



(8) 



(9) 



where e is the energy for the Hamiltonian ([7]). Now expanding (1 + Ax 2 ) 1 for |x| < ^= we can write the 
potential of equation §8§ as 

V(x) = -f +gx 2 - A 0(x 3 ) (10) 
A 

It is clear from (fTO)) that the term (— v) in equation ([9]) cancels from both sides of the equation ([8]), so that the 
new eigenvalues (j9]) are actually the old eigenvalues e of the Hamiltonian (J7]). Also, as A — > 0, the potential 
and the eigenvalues of equation |8]) reduces to those of a linear harmonic oscillator. 

Now generalizing the potential of the equation ([8|) as below, the corresponding PDMSE now reads 



(1 + Ax 2 ) 



d 2 ip 
dx 2 



Xx 



dip 
dx 



B 2 



A 2 -A^X 



1 + Ax 2 



B{2A + VX) 



VXa 



1 + Ax 2 



■ A 



ip = Eip 



(11) 



It is seen from (| 1 1 1) that if we put B = then the potential reduces to that of the nonlinear oscillator with 
4 = A 2 + A^/X. It is to be noted that this generalization should correctly reproduce the A — > limit, in which 
case the equation ([TTj) reduces to the Schrodinger equation for linear harmonic oscillator. In Appendix I we 
have shown that in the limit A — > and for A = ^3= (which is one of the solution of the quadratic equation 
A 2 + A\fX = j), B = the potential of equation (fTTj) , the energy eigenvalues (fT5]) and the wavefunction 
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given in (|19|) reduces to those of a linear harmonic oscillator. This particular generalization is made so that it 
corresponds to the hyperbolic Scarf II potential [55] in the constant mass case. In order to solve (jlll) , we now 
perform a transformation involving change of variable given by 



d.v 



— smh 



where 



F(x) = 1 + Xx 2 , A > 
Under the transformation ([1"5|) , Eqn. (TTTj) reduces to a Schrodinger equation 

d 2 tp 



dz 2 



+ V(z)ip(z) =Eip(z) 



where the potential V{z) is given by 

V(z) = (B 2 -A 2 - AV\)sech 2 (zV\\ + B(2A + VA) tank (zVxj sech [zV\ 
The potential (I15[) is a standard solvable potential and the solutions are given by [55] 



A 2 



^ n ( Z )=N n i n (l + sm/l 2 (z\/A)) 2 e - rtan ^sinhizVX)) p(?r s i ir-s-h) L sinh ( zy fi\ 

where N n is the normalization constant , s = r = -M= and Pi a '^\x) is the Jacobi Polynomial 

normalization constants N n , n = 0, 1, 2, ... are given by [23j . 



— -i 1/2 

V A n\ (s — n)T(s — ir — n + 5)T(s + ir — n + ^) 



7T 2- 2s r(2s-n + l) 



The eigenvalues E n are given by 



E n = nV\(2A- nVX), n = 0, 1,2, ... < s 



(12) 
(13) 

(14) 

(15) 

(16) 
The 

(17) 
(18) 



Subsequently by performing the inverse of the transformation (|12p we find the solution of PDMSE (TTTj) as 

1/2 



\/A nl (s — n)T(s — ir — n + ^)T(s + ir — n + |) 
7T 2" 2s r(2s-n + l) 

i n (l + \x 2 )-ie- rtan ~ 1{x ^Pi lr ~ s ~^'~ lr ~ s ~^(ixV\) , 



n = 0,l,2,...<« (= A) 



(19) 



At this point it is natural to ask the following question : Are there other solvable potentials corresponding 
to the mass function m[x) = ( 1+ \ x 'i J ? The answer to this question is in the affirmative. The procedure 
to obtain these potentials is similar and so instead of treating each case separately we have presented the 
potentials and the corresponding solutions in Table 1. The first two and the last two potentials in Table 1. 
are actually the generalizations of the nonlinear oscillator potential. Although the other two potentials in 
the Table are not generalizations of the nonlinear oscillator potential, nevertheless they are exactly solvable 
potentials with the same mass function. 
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2.1 Relation between A-deformed Hermite polynomial and Jacobi polynomial, 
Generating function, Recursion relation 

Here we shall obtain a correspondence between the A-deformed Hermite polynomials [21] and Jacobi polyno- 
mials. We recall that the Hamiltonian for nonlinear oscillator is given by [H] 

„2 



2 dx 2 I 1 + Ax 2 



After introducing adimensional variables (y, A) as was done in ref |21j 

A 



y = \/ax 

the Schrodingcr equation Hip — eip reduces to 

1 / 1 A 2^ 2 1 a d 

The eigenvalues and eigenfunctions for A < are p?T] 



A 



1 + A 



1 + Ay 2 



tp = eip 



ip m (y,A) 



H m (y,A)(l-\A\y 2 ) — 
(m+i)-im 2 A , m = 0,1,2, 



(20) 



(21) 



(22) 



(23) 



where 7i m (y, A) is A-deformed Hermite polynomial whose Rodrigues formula and generating function are given 
in (J27J). For A > 0, 

<My,A) = W m (y,A)(l + Ay 2 )-A 
e m = (m + ±) - ±rn 2 A , m = 0, 1, 2 • • • , N A 

where N\ denotes the greatest integer lower than m\(= jj). On the other hand, putting B = and A = 
in the solution (Tl9|) of Eqn. fTTj) . the eigenfunctions of Eqn. ([2T|) can be written in terms of Jacobi polynomial 

as 

My)^N n (l + Ay 2 y^P^^ h '^ i \iyVA) , n = 0, 1, 2 • • • < (A > 0) (24) 



For A < 0, putting B = 0, A = 



in the wavefunction of the 5th entry of Table 1. and using (|20"|) we obtain 



Mv) = N n {l + Ay 2 )-^pi ' ' i} 



(yV|A|) , n = 0,1,2. 



(A<0) 



(25) 



Comparing Eqns. p^)) and (f^5]) and also Eqns. p3"|) and ([M]) . it is possible to derive a relation between A- 
deformed Hermite polynomial 7i n (y,A) and Jacobi polynomial P n a '^\x) as 



Pi * 4 *>(i„VA) = i(^=) «n(»,A), VA 



(26) 



The Rodrigues formula and the generating function for the A-deformed Hermite polynomial 7i„(y,A) were 
given by [2T] 



n n (y,A) 



1 + Ay 2 



(27) 



(-l)n 2 A+^ 

^(t,y,A) = (l + A(2ty-i 2 ))i 

It was shown [21] that the polynomials obtained from the generating function F(t, y, A) with those obtained 
from Rodrigues formula arc essentially the same and only differ in the values of the global multiplicative 
coefficients. We have observed that if the generating function F(t, y, A) is taken as 



00 1 (--) t n 
(1 + A(2ty t 2 ))^ = £ — ji ^pr-Ti-n (y, A) — 

n=0 V2 Tin 



(28) 
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where (a)„ represents Pochhammer symbol given by (a) n — r ^"' ) then the polynomials obtained from the 
above relation are exactly same with those obtained from Rodrigues formula given in Eqn. (|27p . 
Correspondingly the recursion relations are obtained as 

(A(2n + 1) - 2) [2(1 - nA)yH n (y, A) + (A(2n - 1) - 2)nH n -i(y, A)] = (nA - 2)H n+1 (y,A) (29) 

and 

(A(n - 2) - 2) [2(A(2n - 1) - 2)nH n (y, A) - (A(n - 1) - 2)H' n (y, A)] 

= nA(A(2n-l)-2) [2(A(n- 2) - 2)yH' n _ 1 (y,A) - (n- l)(A(2n- 3) - 2)K-a(tf,A)] 



where 'prime' denotes differentiation with respect to y. For A — > Eqns. ([29| and (|30[) give the recursion 
relations for Hermite polynomial [24) . 



3 Shape invariance approach to supersymmetric PDMSE 

Supersymmetric approach to PDMSE [5] may be discussed either by reducing the PDMSE to constant mass 
Schrodinger equation or start with modified intertwining operators consisting of first order differential opera- 
tors. Here we shall be following the later approach. Thus we consider operators of the form 

A = P x -iW(x), A^ = P x +iW{x), P x = J_ i-i-f) (31) 

VMx) V dx J 

We now consider the supercharges Q, defined by 

•-(::)■ *-(:?) 

The supersymmetric Hamiltonian is then obtained as 

(T :J 

where the component Hamiltonians are given by 

pn, ( 1 d 2 ( ml \ d , W 

H PDM = + t \ +W 2± (34) 

m(x) dx \2m z / dx 
The Hamiltonians H± DM are supersymmetric partners and the potentials are 

yPDM = W 2 {x) ± (35) 

It can be easily seen that the following commutation and anticommutation relations 

Q 2 = Q f2 = [Q,H PDM ] = [Q\H PDM ] = 

{<3,Q + } = {Q ] ,Q} = o 



(36) 



together with Eqn.(|33p complete the standard supersymmetry algebra 22. 25j. For unbroken Supersymmetry 
(SUSY), the ground state of H- has zero energy (E^ = 0) provided the ground state wave function ^ ~\z) 
given by (Atp^ = 0) 



ipQ ^ (x) = N exp 



V m (y) W(y)dy 



(37) 
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is normalizable. In this case it can be shown that, apart from the ground state of H- , the partner Hamiltonians 
H± have identical bound-state spectra. In particular, they satisfy 

4+i = 4 +) ,n = 0,l,2-.. (38) 
The eigenfunctions of H± corresponding to the same eigenvalue are related by 

It may be noted here that the superpotential W(x) and therefore the factorization of the Hamiltonian could be 
generated from the ground state solution of the equation. In a remarkable paper [26) , Gendenshtein explored 
the relationship between SUSY, and solvable potentials. The pair of potentials V±(x,ao), ciq being a set of 
parameters, are called shape invariant if they satisfy the relationship [5J [55] 

V+(x,a ) = W 2 (x,a ) + W'(x,a ) 

= W 2 {x,ai)-W'(x,ai)+R(a ) (40) 
= V-(x,ax) + R(a ) 

where a\ is some function of clq and R(ao) is independent of x. When SUSY is unbroken the energy spectrum 
of any shape invariant potential is given by [22] 

n- 1 

4~ )= I>(ai) ; 4~ )= (41) 

i=0 

We are now going to study the factorization and the shape invariance property of the potentials for the PDMSE 
Schrodinger equation. As an example let us consider the generalized nonlinear oscillator of section [5] For this 
it is now necessary to choose the superpotential W(x) so that i2_ can be identified with the Hamiltonian of 
Eqn. lfTTj) . In this case we choose the superpotential to be 

W = A ^ X +B - 1 (42) 

VI + Xx 2 VTTJx 2 

Therefore the Hamiltonians hP dm and H^ DAI can be factorised as 



H PDM = A \ A 

= -(1 + Xx 2 )£ - Xx£ + B \i\~J^ + B(2A + VA)(^) + A* 

H PDM = AA \ 

= '(I + ^ 2 )£ - Ax £ + *'-4 a + B(2A VX)(^) + A* 



(43) 



These two Hamiltonians are related by 

jj-pdm ( X . A}B ) = h pdm (a;; ^ _ VA, B) + VA (2 A - %/A) (44) 
so that they satisfy shape invariance condition 

Hl DM {x, oo) = H PDM ( Xl a x ) + R(a ) (45) 
where {a } = (A, B), {a x } = (A - V% B) and R(a a ) = \/X(2A - y/X). 
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The ground state ipo(x, ao) of the Hamiltonian H PDM is found by solving Atpo(x, ao) = 0, and has a zero 
energy i.e. 

H PDM ( Xl a )Mx,ao) = (46) 
Now using (|4"5"1) we can see that ipo( x > a i) is an eigenstate of H PDAI with the energy Ei = R(a ), because 

H PDM {x, a )ip {x, ai) = H PDM (x, ai)ifo(x, ai) + i2(ao)Vto(», 01) ^ 
= i?(a )V'o(a;,ai) , [using ([46])] 



Next, using the intertwining relation H PDM (x, ao)A'(x, a ) = A'(x, a )H PDM (x, ao) and equation ([25]), we 
see that 

H PDM (x, a )A^(x, a )Mx, o x ) = ^(or, oojfl^fo aoMfo o x ) = ^ [H PDM (x, ax) + R(a )] ifo(x> a i) 

(48) 

and hence using (146(1 we arrive at 

H PDM { X , a )A^(x, ao)ip (x,ai) = R(ai)A^ (x,a )^o{x,ai) (49) 

This indicates that A>(%, ao)il>o(x, a\) is an eigenstate of H PDM with an energy E\ = R(ao). Now iterating 
this process we will find the sequence of energies for H PDM as 



E { - ] = R i a i) = nV\(2A - n\f\) , E ( Q } = (50) 



i=0 

and corresponding eigenfunctions being 



^„(a;,ao) = A^(x, a )A j <(x, ai).. (x, a n — i )i/}o(x,a n ) (51) 
where a* = /(a,_i) - /(/(.. ..(/(a ))) = f A - *V% and i?(a,) = VX [2 f A - zVa) - VA . 



t times 

We have found a number of other potentials which are shape invariant for the same mass function. For 
all these potentials the energy, wave functions and other parameters related to shape invariance property are 
given in Table 1. 



3.1 Shape invariance approach to PDMSE with broken supersymmetry 

When supersymmetry is broken neither of the wave functions ip^ix) ~ exp[± j x \/ m(y)W(y)dy] are normal- 
izable and in this case all the energy values are degenerate i.e, H + and -ff_ have identical energy eigenvalues 
[221 [26] 

E { - ] = 4 +) (52) 

with ground state energies greater than zero. So far as we know, little attention has been paid till now to study 
problems involving broken SUSY in the case of PDMSE. Broken supersymmctric shape invariant systems in 
the case of constant mass Schrodinger equation has been discussed in ref [57] . Below we illustrate the two step 
procedure discussed in [28] for obtaining the energy spectra in PDMSE when the SUSY is broken. For this, 
we consider the superpotential as 



W(x,A,B) = Ay/\\\ * - \ 3 V1 + XX '\ 0<x<^=, A<0. (53) 
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Then the supersymmetric partner potentials are obtained using |35|) as 

-(A + Bf 

(54) 



The ground state wave function is obtained from (|37|) as 

tP ( ~ ] + (55) 

For A > 0, B > the ground state wave function V>o ' is normalizable which means the SUSY is unbroken. 
But for A > 0, B < and A < 0, £? > , neither of ip^ are normalizable. Hence SUSY is broken in both 
cases. 

We shall discuss the case A > 0, B < 0. In this case the eigenstates of V±(x, A,B) are related by 

ip^(x,a Q ) = A(x,ao)^\x,a Q ) 

ip { ~\x,a ) = Ai(x,ao)ipW{x,ao), (56) 
£(-)(ao)=4 +) M 



Now we can show the potentials in equation (|54|) are shape invariant by two different relations between the 
parameters. 

Step 1 

The potentials of equation (j54|) are shape invariant if we change A — » A + \/jA[ and B ^ B + -^/jXf.The shape 
invariant condition is given by 

U+(x,A#) = U_ (x,A+ ^/\X\,B + v/|A[) + (a + B + 2v/|A[) 2 - (A + B) 2 (57) 



Now for B < 7= it is seen that the superpotential (|53| resulting from change of parameters as above falls 



in the class of broken SUSY problem for which Eq ' ^ 0. Though the potentials of equation ([54]) are shape 
invariant but we are unable to determine the spectra for these potentials because of the absence of zero energy 
ground state. 

Another way of parameterizations A — > A + y/\X\ and B — ► —B gives us 

V + (x,A,B) = V- (x,A+y/\X\,-B) + (a-B + ^\X\) 2 ~{A + B) 2 (58) 



which shows that V- and V+ are shape invariant. This change of parameters (A — ► A + y]A] and — * — £>) 
leads to a system with unbroken SUSY since the parameter B changes sign. Hence the ground state energy of 
the potential V- (x, A+y/\X\, —B) is zero. From the relation |58|) we observe that V+(x, A, B) and V- (^A + i/fXl, —B 
differ only by a constant, hence we have 



1>+(x,A,B)=il>-(x,A+y/\\\,-B) 

25<+> (A, B) = (x, A+^\X\,-B) + (A-B + ^\X\) 2 -(A + Bf 



Thus, if we can evaluate the spectrum and energy eigenfunctions of unbroken SUSY H^ DM (x, A + y|A|, —B), 
then we can determine the spectrum and eigenfunctions H? DM (x, A, B) with broken SUSY. In the 2nd step 
we will do this. 
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Step 2 

With the help of shape invariant formalism in case of unbroken SUSY for PDMSE (See section [3]) we obtain 
spectrum and eigenfunctions for V-(x, A + i/]A[, —B) as 

Ei-HA+^l-B) = (A-B+^\+2n^\) 2 -(A-B+^\m) 2 
^\x,A+^/\\\,-B)(xx^(l + \x 2 ) 2 ^iP n Vlxl ' V|A| (l + 2\x 2 ) 
Now using ([60]) , (|59]) and (|56|) we obtain spectrum and eigenfunctions for V~(x,A, B) with broken SUSY as 

£?<->tA,B) - (A- B + ^/\X\ + 2ny/\X\y - (A + B) 2 

1=5. i_ ( L ~ B A - 1 ) (61) 

^-\x,A,B)(xx^{l + Xx 2 ) 2 ^Py 2 J(l + 2\x 2 ) 

Similar approach can be applied in case of A < and B > 0. In this case we change (A, B) into (— A, B + y|Af) 
and the shape invariance condition is 

V+(x,A,B)=V- (x,-A,B + y/\\\) + {B- A+^\\\) 2 -{A + Bf (62) 

And 

E(-\A,B) = (B - A+ ^/\X\ + 2ny/\X\y - (A + B) 2 

B 1 ^ 1^) ( 63 ) 

A, B) oc i (1 + As 2 ) 2 v^ ^ 2j (l + 2Xx 2 ) 

4 Exactly solvable VT symmetric potentials in PDMSE 

Here we shall find exactly solvable complex potentials, some of which are related to the nonlinear oscillator 
potential, within the framework of PDMSE. Before we consider any particular potential let us note that a 
quantum mechanical Hamiltonian H is said to be VT symmetric [2 8) if 

vth = nvr (64) 

where V is the parity operator acting as spatial reflection, and T stands for time reversal, acting as the 
complex conjugation operator. Their action on the position and momentum operators are given by 

V : x — > —x, p — > — p, T : x — > x,p — > —p, i —* —i (65) 

For a constant mass Schrodinger Hamiltonian the condition for VT symmetry reduces to V(x) — V*(—x). 
However in the case of position dependent mass an additional condition is required. To see this we note that 
in the present case the Hamiltonian is of the form 

^ Id 2 m'(x) d ^ 

2m{x) dx 2 2m 2 (x) dx 

From (|6"5)) it follows that the conditions for the Hamiltonian (ffi6"|) to be VT symmetric are 

m(x) = m(-x) , V{x) = V*{-x) (67) 

It may be pointed out that here we are working with a mass profile m(x) — (1+Ax 2 ) -1 which is an even function 
and consequently satisfies the first condition of (|67[) . To generate non-Hermitian interaction in the present 
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case we introduce complex coupling constant. As an example let us first consider the potential appearing in 
(fTTj) . It can be seen from (fT5|) that the energy for this potential does not depend on one of the potential 
parameters, namely B. Thus we consider the complex potential 



V{x) = 



B 



a 2 -aVx 



1 + Xx 2 



iB(2A+Vx) 



l + Xx 2 



A 



(68) 



From ([68]) it can be easily verified that V(x) = V*(— x) so that the Hamiltonian (|66| with this potential is VT 
symmetric. In this case the spectrum is real and given by l|18[) . Proceeding in a similar way we have obtained 
the spectrum of a number of VT symmetric potentials and the results are given in Table 2. Incidentally all 
the potentials in Table 2 are shape invariant and the results can also be obtained algebraically. 



5 Quasi exactly solvable PT symmetric potentials in PDMSE 

The complex sextic potential in the constant mass Schrodinger equation has been discussed in ref [29] . By 
using the transformations (112|) for A > we obtain the corresponding quasi exactly solvable potentials in 
PDMSE. 

For A > 0, the potential is taken as 

6 _ k 

V{x) = Y.^(sinh{xVX)) (69) 

where for V(x) to be VT symmetric, c\, c 3 , C5 are purely imaginary and c 2 , C4, cq are real. 
Following ref [25], the ansatz for the wave function is taken as 

i/,(x) = f(x) exp ( - \(sinh(xy/\))- j ] (70) 

where f(x) is some polynomial function of x. We shall focus on the following choices of f(x): 

(a) f(x) = l 



(^sinh(x^/\)j 



(b) f(x) = ± — I +a a 

(sinh{x\/~\)) isinhlx^TX)] 
(c) f(x) = ± J +a a 

For complex potentials, ao is purely imaginary in (b), but in (c) a\ is purely imaginary, but ao is real. 
Without going into the details of calculation, which are quite straightforward, let us summarize our results. 

Case 1: f(x) = 1 

In this case the relation between the parameters Cj and bi are found to be 
ci = -3b 3 + 2b 1 b 2 , c 2 = -6b 4 + 3b 1 b 3 + 2b 2 , c 3 = 4b 1 b 4 + 6b 2 b 3 c 4 = 8b 2 b 4 + h 2 3 , c 5 = 12b 3 b 4 , c 6 = 8b 2 4 (71) 



and 



E = b 2 - \b\ (72) 
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Without loss of generality, we can choose cq = \ which fixes the leading coefficient of V[x). It gives 6 4 = ±i. 
Taking the positive sign to ensure the normalizability of the wave function we obtain 



ip{x) = exp 



( 



\ 



b\ (sinh(xVX)j 62 ( sinh(xy/X)\ 63 (sinh(xy/X)\ (sinh(x^/X)j 



VA A Aa/A 4A 2 



(73) 



Now if 61 and 63 are purely imaginary then c±, C3, C5 are also purely imaginary. In that case V(x) in Eqn. 
and ip(x) in Eqn. ([70]) are VT symmetric and E is real. 

sin 

Case 2: f(x) — — -= — \- ao, where ao is purely imaginary. 

yA 

In this case wave function is of the form 



exp 



(sinh(x\/~X) 

bi (^sinh(x\f\)j 62 (^sinh(x\/\)j 63 ^sinh(x\/\)j (sinh(x-\/X) 



ip(x) = I -= h a 



(74) 



1 



V\ A AVA 4A 2 

In this case the relation between the parameters are given by 

5 9 
ci = -6fe 3 + 26i6 2 + ao , c 2 = -- + 36163 + 26 2 , c 3 = 61 + 66263C4 = 26 2 + -6g , c 5 = 36 3 , c 6 = | (75) 

ao satisfies the condition 

a 3 ~ 3b 3 a 2 + 26 2 a - 61 = (76) 

The energy is given by 

E = --b\ + 36 2 - 3a & 3 + ao ( 77 ) 

We now consider two special cases, 
(a) 61 = 63 = and a 2 , < 0. 

In this case c\ is purely imaginary and c 3 = C5 = 0. Moreover c\ = ao = ±i\/262- So we get two different 
complex potentials corresponding to above two values of c\ with same real energy eigenvalues. The potential, 
energyvalues and the eigenfunctions are given by 

V( X ) = 1 (sinh(xV\))-° + ifr( ainh ( xy /xfi-4. + <&jJH(8inh(xy/\))- 2 ± iffi (sinhixVX))' 1 
E = 6 2 > (78) 
= ^mMjVl)- 1 ± l ^2b^j exp (mnh(xVX)) ~ % - ^ (sinh{xV\)) 

It can be easily seen from the above equations that the potential is VT symmetric, while the wave function is 
odd under VT symmetry. 

(b) 61 =0,63 ^0 
Then from ([75]) we get 



^ = -(363 +^961- 862) (79) 
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So in order to make oo imaginary we must have % 2 — 862 < 0. or b 2 = — | £>3 1 2 < |6 2 . 

In this case also there exist two different complex potentials corresponding to two values of 63 with the same 
real energy eigenvalues E = 3&2 — 3ao&3 + a 2 ,. 



Case 3: f(x) = (™^)" 



a 1 



7t 



do, where a\ is imaginary and ao is real. 



In this case the relation between the parameters is given by 



a x = 2b 3 



a 



2bo 



b z 3 ±J(2b 2 -3biy + 2 



(80) 



The wave function, energy and the potential are of the form 



(sinfrCsVA))- 2 + ^ (.inM^A))- 1 + 1 ^ _ fe 2 ± ^ _ 352)2 + 2) 



exp 



-2b 3 (b 2 - b%) 



(sinh(x^/X)\ ^sinh(x\/X)j (^sinh(x\f\)^j ^ fswi/i(:E\/A)l 

7t 62 A 63 AVI 4 A 1 



£± = -2^3(62 - bf) 2 + 3b 2 - &§ ± J(26 2 - 3& 2 ) 2 + 2 



(81) 
(82) 



V(sb) = -^(szn/i(a;\/A))- 6 + -^=(«n/i(a:VA))- 8 + 



3& 3 



(26 2 



2b 



2A 3V v " A 2 \/A 
[2(^+3&2^-3&|)-I] 

A 



A 2 



^(smft^))" 4 + ^L(4& 2 - ^)(szn/i(xV / A))" 



A\/A 



(sinh(xVX))- 2 + HAh2 ~ A ^ bl ~ 7 \ sinh{xJ\))- 1 

v A 



(83) 

The result (|80|) to (|83|) are valid both for real and purely imaginary bi. When bi are purely imaginary the 
potential and wave function are VT symmetric while for real bi VT- symmetry is broken. In particular when 
&3 is purely imaginary we have a complex VT- symmetric two parameter family of potentials corresponding 
to two values of ao with two distinct real eigenvalues. 



6 Discussion 

We have studied various exactly solvable as well as quasi exactly solvable and non-Hermitian generalizations 
of the quantum nonlinear oscillator with the mass function ( x+\x 2 ) • We nave a ^ so obtained a closed form 
normalisation constant for the eigenfunctions of quantum nonlinear oscillator. A relationship between the A 
deformed Hermite polynomial and Jacobi polynomial has also been found. By exploiting supersymmetry of 
the PDMSE we have obtained some shape invariant potentials corresponding to this particular mass function. 
We have considered the shape invariance approach to PDMSE with broken supersymmetry as well. As for 
the future work we feel it would be interesting to examine Lie algebraic symmetry of the exactly solvable 
potentials. In view of the fact that in the present case the transformation (|12[) is invertible, it seems promising 
to study whether or not the Lie algebraic symmetry of the constant mass system can be transported back to 
the non constant mass case. Another interesting area of investigation would be to study the classical analogs 
of some of the models (especially the VT symmetric ones) considered here. 
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Appendix I 

For B = 0, A = the potential of equation (fTTj) and it's energy eigenvalues (fl8|) reduces to 



V(x)= ^---ajil + Xx 2 )-^^- (II) 

E n = 2na - nX (12) 
For |ac| < the potential (II) can be written as 

a 2 

(13) 



2 



V (x) = ^~ - a) (1 - Ax 2 + AV - A 3 x 6 + ...) + j- 

= a 2 x 2 — \{a 2 x A — Aa 2 s 6 + ...) + \{ax 2 — Xax 4 + ...) — a 

For A — > the potential reduces to 



V(x) = a 2 x 2 - a (14) 

It is clear from (14) and (12) that for A — * the potential (JTTJ) and the energy eigenvalues (TT5)) reduces to 
those of a simple harmonic oscillator. 

For A = ,B = Q and using the relation (|2"6")l the expression for the wavefunction (fT^|) is as 



V„(z) = N' n (l + Xx 2 )-^ H n (y^x, ~) 



(15) 



where 

1 /ax § 
<=2^(a) 



""(f-") r (f~"+T) r (f-"+^) 
tt«! ! ! " _ f A""ir(^-n+l) 



1/2 (16) 



Now for A — > the A-deformed Hcrmite polynomial becomes the conventional Hermite polynomial H n [21]. 
Consequently at A — ► limit the unnormalized wave function given in equation (15) reduces to 

ax 2 

ip n (x) tx e ~H n (y/ax) (17) 
Using the asymptotic formula T(az + b) ~ \/2Tre~ az (az) az+b ~2~ (see 6.1.39 of the ref[24]) in (16) we have 

(/— TtA \ 1 / 2 
( I8 ) 

Therefore from equations (17) and (18) it follows that for A — ► the wave function given in equation (fT9|) 
reduce to that of simple harmonic oscillator. 
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3 2 -A 2 -AVA +b(2 A+ yT)^Sa_ + A 2 



- h 



ix/A(2A - rav'X) 



'(l + A^) 2 



,-r tan- 1 (xVX) 



(A - ./I, B) 



VT I2A - (2t + l)vTl 



; («*Va) 



4 2 + 4- A( - 4 + V ?' ) + 2B^2i- , B < A 2 

A? ^/l + Ax 2 



i -^Aa + . 

\/l + Ax 2 



• 42 + 45- - f- 4 - "^) 2 



\/l + Ax 2 / I \A + Ax 2 



(A - B) 



,("l. s 2) I ivT 



\/l + Ax 



^ " 2 ( yTTH 



a 2 - [a - (i + i)vT] 2 

+ B 2 _ B 2 

[A-(,+ l)v/T] 2 



VAa: Ax' 1 
< ^VA" < oo 



(A + n V\) 2 



■ + A 



(A + iv^X, B) 



P] («3» s 4) 



A 2 +B 2 + Av^ _ ~ " x) £ l + >;f +A 2 ~ 



/ l + Ax 2 A 



A 2 - [A + (i + 1)VA] 2 

+ B 2 B 2 

[a+(,+ 1)vT]' 



xV^(2A - nVT) 



< xv'A < 



1 + Ax 2 

(r- 



—4) 



Ax 2 + 1) k 2 



(^) 



(A - ivT, B) 



VT [2 A - (2t + l)vT] 



A^+B^-A^A^ 
1 + Ji 2 



- B(2A - v^A|)- 



/|A|(2A + nv^XT) 



(l-xV7A|)' 2 



"'(l + xyiAl)* 2 > 



(a + ivTXT, b) 



/IaT [2A + (2. + ovTXT] 



(*vTaT) 



1 + Ax 2 1+A 



^1 — d 2 _l B 2 



\/l + A^ 



A^ + (A + n^\MV 



(A + n ^lX\) 2 



V 1 + Ax 2 > 



(A + »^|A|, B) 



XT - x " vTxT 



\/l + Ax 



^ [a +( , + i )v ^] 2 
-a 2 + [a + (i + ovTXTl 



Table 1: Exactly solvable shape invariant potentials V(x), superpotential W(x) , energy eigenvalue E n and wave functions ip n {x). Where s = -^,r = 

-j-, n — y, a — si = s — n + a, S2 = s — n — a, S3 = a — n — s, sa = — (s + n + a), s' = —4= & r' = ? .The first four entries correspond to A > and 
vA s n y|A| v|A| 

the last two correspond to A < 0. 



\+L^ A ^ + iB(2A + ^'l+Tx^ + - 



*n(x) 



+ Ax 2 )~f c- r tan I (^V^) 



<J 1 + Ax 2 y' 1 + Ax 



A 2 _ B 2 A(A+VX) 
"A 7 1 + Ax 2 



(A+vOT) VA X 



B < A 2 



\jl+\x 



■ + 



i2 -^-^-"^» 2 +(A^W 



cV\_ 



o(«l.«2) I X vT 

, ^1+Ax 2 



v' 1 + Ax2 , A(A-vT) 

7T 1 ; — ? 5 > 



V^x 
< xVA < 



, B 4 •/ 1 + Ax 



(A + nVI) 2 + 



(A + nVA) 2 



Vl + Ax- 

xvT 



4^- 



p(s 3 ,s 4 ) / 1 + Ax 



x V A x V A 



( 1 + Ax 2 N 

V xvT ) 



A 2 -B 2 +As/\ 



iB(2 A + A)Vi±^ 
< xVA < oo 



■A A , A < B 



nv / A(2A - n-SX) 



(\/l + Ax 2 - 1) ( r 2 S) (^/l + Ax 2 + 1) <r 2 S) 
n(r _ s _l,_ r _ s _l ) ^^ / — — 



s' —r' 

(i - *vW) (iL ~^ 1 ~ ) (i + ^VW)^' 
p' s '- r '-i' s ' +r '-^(^) 



a 2 -b 2 -aVW\ _ tB{2A _ yrxp X ^TM _ A 2 



y/l + Ax 2 V 1 + A:l: 



iV^(2A + n^) 



1 < x < 1 



^'^VW' - 2»B "V^T - A 2 - si 
1+Ax 2 1+Ax 2 A 2 



A < x < 

-TXT - - VTM 



\/l + Ax 2 



A 2 + (A + n/W) 2 + 



(A + rev ^) 2 



-( S 2 ") e _ av ^x 



„(-«'— n — <o, — »'-n+io) | , xy^J 

1 yr+r^ 



Table 2: Exactly solvable PT Symmetric potentials, where s = -^,r = , ri = i-f ,a = = s — n + a, S2 = s — n — a, S3 = a — n — s,S4 

— (s + n + a), s' = —7= & r' = i f . The first four entries correspond to A > and the last two correspond to A < 0. 
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